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1. Introduction 

Baryons containing heavy quarks provide an interesting laboratory for studying QCD. They 
combine two different regimes: the slow relative motion of the heavy quarks with the relativistic 
motion of a light quark. The energies, detectors and luminosities at disposal in experiments made 
possible the observation of heavy baryons with one heavy quark. The first heavy baryon signal 
ever seen was the A+ at BNL in 1975 followed soon after (1976) by the discovery of E++ at 
FNAL. Also at Fermilab, in 1981, the first baryon containing a bottom quark A^ was observed. 
Substantial progress has been made, and currently, there is a total of 19 charmed and 5 bottomed 
baryons in the PDG summary tables [1]. In the last decade, charmed baryons have mainly been 
observed at the B -factories, whereas new bottom baryons were found at the Tevatron and more 
recently also at the LHC. Masses, decay form factors, lifetimes and widths have been determined. 
However, identification of the spin and parity quantum numbers from experiments is still missing; 
so far, they are assigned based on quark model expectations. It is expected that the large statistics 
provided by the LHC will allow for their identification through the study of angular distributions 
of the particle decays. Also, the PANDA experiment at the FAIR faciUty and the KEK Super-B 
Factory will look for charmed baryon signals. For doubly heavy baryons, the experimental situation 
is less favourable. In 2003, Selex [2] published evidence for H+ , but no other experiment has so 
far confirmed this channel. 

Inspired by the experimental activity, many theoretical approaches have been used that try to 
reproduce the existing spectra and predict new states: for example quark models [3-11], QCD 
sum rules [12-15], heavy quark effective theory (HQET) based models [16], and lattice QCD 
(LQCD) [17-24]. 

In this work, we focus on the study of singly and doubly charmed baryon low lying spectra 
including positive and negative parity states using 2+1 light dynamical flavours. This is achieved 
through the Monte Carlo evaluation of a path integral in discretised Euclidean space time with a 
lattice spacing a. Systematic errors are under control if the ultraviolet cutoff, a~^, is larger than 
the scales of the problem, the size of the box, L, is larger than the typical size of the hadrons under 
consideration and an extrapolation to the physical quark masses is performed. Considering the 
hierarchy of the quark masses and energy scales, for current lattice spacings, 

miight < Ms ^ A < rric < l/a < Mb 

where A is the typical hadronic scale. The charm quark mass is below the cutoff and a relativistic 
calculation of the baryon spectra containing charm quarks is viable. However, it is important to 
show discretisation errors are under control. 

This write-up is organised as follows. In section 2, details of the computational setup are given. 
In section 3 we briefly describe the methodology used to extract the energy levels. Following this, 
a detailed list of interpolating operators is given and results are presented in section 4. We finish 
with some concluding remarks in section 5. 

2. Computational details 

In our calculations, two different sets of gauge configurations with u/d and s sea quarks were 
used, 2-HEX [25] and SLiNC [26, 27]. They were generated by the BMW-c and QCDSF coUab- 
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orations, respectively. They were both generated using a gauge action where 0{a^) effects were 
reduced. The 2-HEX configurations employ tree-level clover Wilson fermions coupled to links 
with two steps of HEX smearing [28]. This means that the lattice artifacts are reduced to O(asa). 
SLiNC configurations use non-perturbatively clover Wilson fermions with stout link smearing on 
the derivative terms [29]: 0{a) effects are non-perturbatively removed. 

The 2-HEX ensembles were generated for a range of lattice spacings from a ~ 0.092 fm down 
to a ~ 0.054 fm and pion masses from M„ ~ 520 MeV down to Mj^ ~ 120 MeV. The spatial 
dimensions cover the range 1.7 — 5.9 fm and the number of configurations per ensemble used for 
this study is ~ 200. See [25] for more details. With this set of ensembles we can perform a 
controlled chiral and continuum extrapolation. 

For the SLiNC configurations, there is only one lattice spacing available, a ~ 0.0795 fm. The 
fight quark masses were tuned to the SU(3)flavour-symmetric point, where the flavour singlet mass 
average niq = (m„ + mj + mv)/3 takes its physical value. Then, and are varied keeping 
niq constant [26, 27]. This is motivated by the Gell-Mann - Oakes - Renner relation as well as by 
SU(3) chiral perturbation theory (;|^-PT): mk approaches its physical value from below, so that it 
is within the range of applicability of X-VY {niK < 600 MeV). Ensembles are available for two 
spatial extents, 1 .9 fm and 2.5 fm. So far, we have analysed one of the smaller configurations, with 
M„ = 348 MeV and Mk = 483 MeV. 

3. IMethodology 

Let ^1, ^2 be two interpolating operators that overlap with the baryonic state we are interested 
in. The correlation function of S'l and S'2 is given by. 



Cp=o(^i,4,0 = (^1(0)4(0) =^lini^Tr 



'£mi\n){n\A\0)e-^"', Z{T) =Tr 



(3.1) 



where T is the temporal extent of the lattice. Every correlation function contains contributions 
from all states with the same quantum numbers. Since we are interested in the low lying spectra, it 
is desirable to eliminate the contribution of excited states. This can be done using the variational 
method [30,31]. 

3.1 Variational method 

The idea is to choose a basis of operators ffi with different overlaps with the state we are 
interested in. A suitable possibility is to create operators with the same Fock structure and different 
spatial extent. This can be achieved by applying different numbers of steps of Wuppertal smearing 
[32,33] to the fermionic fields. A cross correlation matrix can then be built, [C(?)]j-y = C{S'i, S'j,t). 
Solving the generafised eigenvalue problem (GEVP), 

C-i/2(^o)c(OC-i/2(^o)v«(?,?o) = A«(?,?o)v«(?,?o). (3.2) 

It can be shown that the eigenvalues have the following behaviour, 

A«(?,?o)oce-'E«('-fo) [i+o(e-^«')] , A£« with a' > a. (3.3) 
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In this way, the lower levels can be extracted cleanly. Nevertheless, there is some freedom of choice 
of the structure of the interpolating operators. We now motivate the choices made in this work. 

3.2 Interpolating operators 

To construct sensible interpolating operators overlapping with the baryon states, it is useful 
to think in terms of their possible inner structure. The constituent quark model guided by the 
approximate SU(2) or SU(3) flavour symmetry is quite successful when applied to the light quark 
sector. This statement may no longer be true when we consider heavy baryons. As far as singly 
charmed baryons are concerned, predictions of HQET tell us that the light degrees of freedom form 
a diquark and move around the approximately static heavy colour source. In the limit mq — )• oo, the 
total angular momentum of the light diquark becomes a good quantum number. Its total spin 
takes two possible values, and 1 corresponding to a flavour antisymmetric or symmetric structure, 
respectively. It is possible to construct operators following these prescription, cf . Table 1 . 



Singly charmed 



J = 



(S) (I) 


Sd 


(qq)Q 




Name 


Name 


(0) (0) 


(0)+ 


(n')c 




= eahc{l"^Cr5l'''y 


A, 




(0) (1) 


(1)+ 


{ll)c 






Z, 


y * 


(-1) (i) 


(0)+ 


{ls)c 




= E^i^ii^^crss'^y 






(-1) ik) 


(1)+ 


{ls)c 






'~-'c 




(-2) (0) 


(1)+ 


{ss)c 








a* 



Table 1: Summary of quantum numbers of singly charmed baryons and interpolating operators following the HQET 
approach. / and /' stand for light quarks and 5,/ and i-^ stand for strangeness, isospin, and diquark total spin, respectively. 
These operators were suggested in [17] for lattice calculations. 

As we can see, operators (with = 1) have contributions from states with total spin ^ and 
|. We thus have to use spin projectors to disentangle them. For zero momentum, these projections 
amount to 



33/2 ^ 



8'^ 
1 



1 



/y, /,7 = {1,2,3}, 



(3.4) 



Since flavour symmetry is not a good symmetry for charmed baryon systems, in principle it 
is not clear that interpolating operators falling into the irreducible representations of SU(4)flavour 
have a decent overlap with the baryon states. However, these states have also been constructed, 
and we have seen, cf. section 4, that they lead to results that are compatible with the energy levels 
obtained from the HQET based interpolating operators. The tensor product of three fundamental 
representations of SU(4)flavour is, 



4 



20s 



20m O 20] 



'M 



We can construct interpolating operators falling into these multiplets. Since there are at most three 
different flavours in a baryon, the structure is exactly the same as for the interpolating operators 
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in the SU(3)flavour multiplets (with u,d,s quarks). We list the operators corresponding to baryons 
containing {ll''">c) as constituents. 

• SU(4) 20-PLET CONTAINING SU(3) OCTETS — ^ 

ILc- = [c"{xY{Cys)l''{x)\ Fy{x). 

Ac : €f^{x) = ^e* {2 (x)^ {Crs^'Hx)] + [c"{xf{Crs)l'''{x)] ^{x) - [c'- [xf (Cy,)!'' {x)] l'^^{x)] . 
E?: 0^{x) = ^E""^- { [l^{x)T{CY5)c''{x)] {^{x) + [i-(xf{Cys)c\x)\ l^,{x)] . 

• SU(4) 20-PLET CONTAINING SU(3) DECUPLET | | 

Z* : fff' = {2 (c"^(Cr^)/'') /^+ (/"^(C7p)/'') 4} . 

The operators for the baryons with strange content are obtained by substituting the / quark by 
s in the operators given above. 




Figure 1: Structure of a doubly heavy baryon: HQET picture (left), Quarkonium-like picture (right), 

Next, we discuss possible interpolating operators for the doubly charmed baryons. In Figure 
1, we depict two possible structures for these systems. On the left hand side, the HQET picture, 
the diquark is formed by the two heavy quarks QQ and interacts with the remaining light quark as 
if it was a heavy light meson. In this case, the radius of the QQ system is much smaller than A^^ 
The operators reflecting this structure are listed in Table 2. 



Doubly charmed = 3 = i 



(S) (I) Sd 


(QQ)q 


e 


Name 


Name 


(0) (0) (1) 


{cc)l 








(-1) {\) (1) 


{cc)s 









Table 2: Summary of quantum numbers of doubly heavy baryons and interpolating operators following the HQET 
approach. 1,1' ,S, I and are the same as in Table 1. 

On the right hand side, the quarkonium-like picture, the diquark is formed by a heavy and a 
Ught quark (Qq) and will interact with the remaining Q as if it was a Q. In this case the radius of 
the QQ system is larger than ' . An example of these operators is given by 

Once the correlation functions are constructed out of two interpolating operators, one has to 
project the correlator into the desired parity, i.e., 

C{^i, ^J,t) = Tyj{ffi,j{0)J~y{t)), (3.5) 

where T is a polarisation matrix that projects onto positive or negative parity. 
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4. Results 

4.1 Continuum limit 

In Figure 2, we present preliminary results for the continuum limit extrapolation of the Q.* and 
the D.CC using results from the 2-HEX configurations. A combined fit including the lattice spacing 
and the light quark mass dependences has been carried out. The following fit function to the baryon 
masses has been used, 

/iT(a,M;„M,,|A) =/™7x(a) (1 +A2X„) (1 +Asx,) , (4.1) 
where Mss refers to the strange-strange pseudoscalar and 



For the continuum extrapolation, two fit functions have been used, 

fi{a) = 1 +Aia^, f2{a) = 1 +A2aa,. 



(4.2) 



(4.3) 



/i is also considered since one expects one loop corrections to be small due to the use of smeared 
links in the action. The errors were estimated using the histogram method combined with the 
bootstrap statistical analysis [25, 34, 35]. Different fit ranges were chosen and only those with 
;^^/dof < 2 were included. 

In Figure 2 the results for the masses of (left) and D.cc (right) ai^e shown. To illustrate 
the quality of the fit, the data points are shifted to the physical light quark mass values and then 
averaged for each lattice spacing. 
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Figure 2: Preliminary. Continuum limit extrapolation from the 2-HEX configurations for Q.* (left) and flee (right) 
states. The Triangle (violet) stands for the experimental value. The diamonds (black) correspond to the continuum 
extrapolation and the circles (red) are the simulation data at finite lattice spacing shifted to the physical masses of the 
light quarks. Two fitting functions have been used, O(a^) and 0(o!ja). 



4.2 Spectra 

The charmed baryon spectra for the SLiNC ensemble are shown in Figure 3 (left). A 3 x 3 
correlation matrix has been constructed per interpolating operator, with three levels of smearing. 
We show results for the low lying singly (above) and doubly (below) charmed baryon spectra, 
including negative parity states, compai^ed with the experimental results. Simulations are still at 
an early stage and only one combination of light quark masses has been analysed. We can see 
that the mass differences between baryons containing u,d quarks and the ones with s quarks are 
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smaller than the experimental values. This is not surprising since we are far from the physical quark 
masses: the singlet quark mass niq is tuned to the physical value which means that the u,d and s 
quark masses are heavier and lighter than their physical values respectively, in spite of the fact that 
not all systematics are accounted for at present. On the right hand side, we can see a summary of 
lattice results for the singly (above) and doubly (below) charmed spectra, with different systematics 
in each case [18-24]. Also, the continuum extrapolated points from the 2-HEX configurations are 
included. We can see that, overall, lattice results agree with experiment. 



EXP pos P 
EXP neg P 



GT pos P ' 
GT neg P ' 



HQET pos P > 
HQET neg P > 



QUARKS: (uuc) 
SPIN: 1/2 3/ 



GT pos P <■ 
GT neg P >■ 



(udc) 



HQET pos P I » I 
HQET neg P —m-^ 



♦ i 



QUARKS: " 
SPIN: 1/2 




QUARKS: (ncc) 

SPIN: 1/2 3/2 



Figure 3: Singly charmed (top) and doubly charmed (bottom) low lying spectrum. On the left hand side, results from 
the SLiNC configurations are shown (M^j = 348MeV). The errors are statistical only. On the right hand side, a summary 
of lattice results is presented, including results from this work (2-HEX, black squares). 

5. Conclusions and outlook 

Heavy baryons are good systems to probe QCD dynamics. The last decade witnessed a huge 
experimental progress in the discovery of new singly heavy baryons. In the near future, spin and 
parity quantum number identification will be possible, thanks to the large statistics and advanced 
detectors at the LHC. In the longer term, further progress is expected from planned experiments 
(PANDA at the FAIR facility and KEK Super-B factory). 

In this write-up we have presented preliminary results of an on-going project to obtain the low 
lying spectra of singly and doubly charmed baryons on the lattice, including states with positive 
and negative parity. Two different configuration ensembles are being used, employing SLiNC and 
2-HEX fermions. 

In the near future, we expect to analyse more sets of SLiNC configurations at different quark 
masses and for larger volumes. For the 2-HEX ensembles we will focus on expanding the analysis 
to other states. 
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